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ABSTRACT 

A mathematical approach to investigate particle acceleration at shock waves 
moving at arbitrar y speed in a medium with arbitra ry scattering properties was 



first discussed in (jVietri 20031 ; iBlasi fc Vietri 20051 ). We use this method and 



somewhat extend it in order to include the effect of a large scale magnetic field 
in the upstream plasma, with arbitrary orientation with respect to the direction 
of motion of the shock. We also use this approach to investigate the effects of 
anisotropic scattering on spectra and anisotropies of the distribution function of 
the accelerated particles. 

Subject headings: cosmic rays - shock waves 



1. Introduction 

The theory of particle acceleration at shock fronts moving with arbitrary speeds (from 



newtonian to ultra-relativistic) can be formulated in a simple and exact form (jVietri 20031 ; 
[l|), at least in the so-called test particle limit, which neglects the dynamical reaction of 
accelerated particles on the shock. In this framework, all the basic physical ingredients can 
be taken into account in an exact way, with special reference to the type of scattering that is 
responsible for the particles to return to the shock front from the upstream and downstream 
plasmas. The information about scattering is introduced in the problem through a function 
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w(n, fi') which expresses the probability per unit length that a particle moving in the direction 
y! is scattered to a new direction fi. It is worth stressing that w can have a different functional 
form in the upstream and downstream plasmas, in particular in the case of relativistic shocks. 

The repeated scatterings of the particles lead eventually to return to the shock front, 
as described in terms of the conditional probability P u (fJ>o, aO (Pd(yo, /•*)) that a particle 
entering the upstream (downstream) plasma in the direction /i , returns to the shock and 
crosses it in the direction of the downstream (upstream) plasma in the direction fi. The 
mathematical method adopted to calculate the two very important functions P u and 
based upon the know ledge of the elementary scattering function w was described in detail in 
( IBlasi &: Vietri 20051 ). and is based on solving two non-linear integral-differential equations 
in the two independent coordinates fi and fi. 



(jVietri 20031 ) showed on very general grounds that the spectrum of accelerated particles 
is a power law for all momenta exceeding the injection momentum. The slope of such power 
law and the anisotropy pattern of the accelerated particles near the shock front are fully 
determined by the conditional probabilities Pd and P u and by the equation of state of the 



through different methods, both semi-analytical (see for example ( 


Kirk & Schneider 1987; 


Gallant & Achterbere 1999 


Kirk et al. 2000: Achterbere et al. 2001 


)) and numerical, by us- 


ine Monte Carlo simulations (e.a. (Bednard & Ostrowski 1998; iLemoine & Pelletier 2003; 


Niemiec & Ostrowski 2004; 


Lemoine & Revenu 2005)1. The theorv 


of particle acceleration 
d versus several of these 


developed bv (1 Vietri 2003) and (Blasi & Vietri 2005) has been checke 



calculations existing in the literature, both in the case of non relativistic shocks and for rel- 
ativistic shocks, and assuming small as well as large pitch angle isotropic scattering (see 



( IBlasi fc Vietri 20051 ) for an extensive discussion of these results). 



In this paper we extend the application of this new theoretical framework to two new 
interesting situations: 1) presence of a coherent large scale magnetic field in the upstream 
fluid; 2) anisotropic scattering. In both cases we calculate the spectrum of accelerated 
particles and the distribution in pitch angle (upstream and downstream) for shock fronts 
moving with arbitrary v elocity. The results of point 1) are compared with those obtained in 
( lAchterberg et al. 20011 ). carried out for a parallel ultra-relativistic shock. 



The paper has been inspired by the need to address several points of phenomenological 
relevance. As far as relativistic shocks are concerned, it was understood that the return 
of the particles to the shock surface from the upstream region can be warranted even in 
the absence of scat tering, provided the ba ckground magnetic field is at an angle with the 
shock normal (e.g. (lAchterberg et al. 200ll )). This is due to the fact that the shock and the 
accelerated particles remain spatially close and regular deflection takes place before particles 
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can experience the complex, possibly turbulent structure of the upstream magnetic field. 
This implies that the calculation of the spectrum of the accelerated particles cannot be 
calculated using a formalism based on the assumption of pitch angle diffusion, as in the vast 
majority of the existing literature. 

In the downstream region, the motion of the shock is always quasi-newtonian, even 
when the shock moves at ultra-relativistic speeds. This implies that the propagation of the 
particles is generally well described by (small or large) pitch angle scattering. However, the 
turbulent structure of the magnetic field, responsible for the scattering, is likely to have an 
anisotropic structure and be therefore responsible for anisotropic scattering. In fact, even 
in the case of isotropic turbulence, the scattering can determine an anisotropic pattern of 
particle scattering. It follows that a determination of the spectrum able to take into account 
these potentially important situations is very important. 

The outline of the p aper is the follow ing: in jj2] we briefly summarize the theoretical 



framework introduced in (IVietri 20031 ) and (IBlasi fc Vietri 20051 ). In £j3]we consider in detail 



the case of a large scale magnetic field in the upstream frame and no scattering of the 
particles. The scattering is assumed to be isotropic in the downstream plasma. In §H we 
introduce the possibility of anisotropic scattering in both upstream and downstream plasmas. 
We summarize in 351 



2. An exact solution for the accelerated particles in arbitrary conditions: a 

summary 

In this section we summariz e the main characteristics of the theory of particle acceler 



ation developed by (jVietri 20031 ) and (IBlasi fc Vietri 20051 ). The reader is referred to this 



previous work for further details. The power of this novel approach is in its generality: it 
provides an exact solution for the spectrum of the accelerated particles and at the same time 
the distribution in pitch angle that the particles acquire due to scattering in the upstream 
and downstream fluids. This mathematical approach is applicable without restrictions on 
the velocity of the fluid speeds (from newtonian to ultra-relativistic) and irrespective of the 
scattering properties of the background plasmas (small as well as large angle scattering, 
isotropic or anisotropic scattering). The only condition which is necessary for the theory to 
work is common to most if not all other semi-analytical approaches existing in the literature, 
namely that the acceleration must take place in the test particles regime: no dynamical 
reaction is currently introduced in the calculations. As a consequence, the shock is assumed 
to conserve its strength during the acceleration time, and the acceleration is assumed to have 
reached a stationary regime. 
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The directions of motion of the particles in the downstream and upstream frames are 
identified through the cosine of their pitch angles, all evaluated in the comoving frames of 
the fluids that they refer to. The direction of motion of the shock, identified as the z axis, is 
assumed to be oriented from upstream to downstream, following the direction of motion of 
the fluid in the shock frame (/i = 1 corresponds to particles moving toward the downstream 
section) . 



The transport equation for the particle distribution function g, as obtained in (IVietri 20031 ) 
in a relativistically covariant derivation reads 



in which both scattering and regular deflection in a large scale magnetic field are taken into 
account. 

Here all quantities are written in the fluid frame, with the exception of the spatial 
coordinate z, the distance from the shock along the shock normal, which is measured in the 
shock frame, u and 7 are, respectively, the velocity and the Lorentz factor of the fluid with 
respect to the shock. 9 and are the polar coordinates of particles in momentum space, 
measured with respect to the shock normal, while <fi is the longitudinal angle around the 
magnetic field direction. As usual \i = cos 8 and u = eB/E is the particle Larmor frequency. 
W(p,, fi',(j),4>') is the scattering probability per unit length, namely the probability that a 
particle moving in the direction (//, 0') is scattered to a direction (/i, 0) after travelling a 
unit length. 

An important simplification of eq. ([T]) occurs when an axial symmetry is assumed. In 
this case the scattering probability depends only on A = <fi — <p' and the large scale magnetic 
field can be either zero or different from zero but parallel to the shock normal. In both cases it 
is straightforward to integrate eq. ([1]) over <fi: the two-dimensional integral on the right-hand 
side simplifies to an integral in one dimension, while the term u{dg/d<p) disappears. 



These simplifications lead to 

dz 



dg 

j(u + n)—= I [-w(fj!,fj)g(fj) + w(fj,,fj,')g(iJ, , )]dfj,', (2) 



where 



u>(/i, //) = J W((jl, //, A) dA and 

#0) = 2~ / g{n,<P)d<P- 
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The physical ingredients are all contained in the two conditional probabilities P u (fi ,fi) 
and Pd(fi 0) fi): these two functions provide respectively the probability that a particle en- 
tering the upstream (downstream) plasma along a direction ft exits it along a direction 
fi. In the absence of large scale coherent magnetic fields, the two functions P u (fi , fi) and 
Pditto, u) were defin ed through a set of two integral-differential non linear equations by 
(IBlasi fc Vietri 20051 ). We report these equations here for completeness: 



-P«(/io,/i) 
r—u 

' dfi 



d(fi ) d(/j) \ _ w(fi,fi Q ) 

u + fl J U + fl 



U + fl 
,w(fl' ,fl )P u (fl', fl) 
U + fl Q 



dfi' P u (fi' , fi) / dfi 



-i 



,w(fi,fi)P u (fi ,fi) 
du ■ — + 

U + fl' 



U + fl" 



(3) 



dfi 



d(fi ) d(fi) \ w(fi,fi ) 



u + Ho u + h 

,P d (// o,//)w(//,//) ^ 
U + fl' 



In the equations above we used: 



U + fl a 



+ 



1 ^^(/AaOM/A^o 



du! P d (fi' ', fi) / dfi 



—u 
—u 



U + fl 
„w(fl',fl")P d (fl ,fl") 



U + fl" 



(4) 



+i 



d(fi) = / w(fi', fi)dfi' 



(5) 



which is unity by definition. 



It is worth stressing that Eq. (J3j) provides automatically the correct normalization for 
the return probability from upstream: f_ dfi'P u (fi ,fi') = 1, independent of the entrance 
angle fi Q . In £j3]we will generalize the method to include the possibility of deflection by large 
scale magnetic fields, which is one of the achievements of this work. In that case we will 
show that the return probability from upstream is no longer bound to be unity, due to the 
escape of particles from the upstream region. 

The pro cedure for th e cal culation of the slope of the spectrum of accelerated particles, 
as found by (IVietri 20031 ) and (IBlasi fc Vietri 20051 ). is as follows: for a given Lorentz factor 
of the shock (7,,), the velocity of the upstream fluid u = f3 s is calculated. The velocity Ud of 
the downstream fluid is found from the usual jump conditions at the shock and through the 
adoption of an equation of state for the downstream fluid. 

Once th e two functio ns P u and Pd have been calculated, the slope of the spectrum, as 
discussed in ( IVietri 20031 ). is given by the solution of the integral equation: 



- 6- 



where 



Q T (e,/i) 



-1 \l-U re lU 



3-s 



(6) 
(7) 



Here u r , 



rel 



1— UUd 



is the relative velocity between the upstream and downstream fluids 
and g(fi) is the angular part of the distribution function of the accelerated particles, which 
contains all the information about the anisotropy. Note that in Eq. ([7]) all variables and 
functions are evaluated in the downstream frame, while the P u calculated through Eq. ([3]) 
is in the frame comoving with the (upstream) fluid. The P u that need to be used in Eq. (J7J) 
is therefore 



djj, 



1 - u 2 rel 



(1 - U re i/J,) 2 _ 

The solution for the slope s of the spectrum is found by solving Eq. ([6]). In general, this 
equation has no solution but for one value of s. Finding this value provides not only the 
slope of the spectrum but also the angular distribution function g(fi). 



2.1. The special case of isotropic scattering 

No assumption has been introduced so far about the scattering processes that determine 
the motion of the particles in the upstream and downstream plasmas, with the exception of 
the axial symmetry of the function W(fi,fi', 



A special case of this symmetric situation is that of isotropic scattering, that takes place 
when the scattering probability W only depends upon the deflection angle 9, related to the 
initial and final directions through 

cos 9 = + y/l — ji 2 a/1 — /i' 2 cos(0 — 0') . (8) 

Among the many functional forms that correspond physically to isotropic scattering, 
the simplest one is 

i f 4 „ 4 1 1 — cos © 

W(n,/j',<j),(f)') = W(cose) = -e — , (9) 

a 

where a is the mean scattering angle. Integration of eq. (Q over <fi — <p' leads to 



1 _iV. / v/i-z/V 1 -// 

-e ° In I 



«;(//,//) = -e - J , (10) 

a \ a 
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with Iq(x) the Bessel function of order 0. Eq. (Q, first introduced in (IBlasi fc Vietri 20051 ). 
naturally satisfies the requirement of being symmetric under rotations around the normal 
to the shock surface. In the limit u < 1 this function becomes a Dirac Delta function, 
strongly peaked around the forward direction, corresponding to isotropic Small Pitch Angle 
Scattering (SPAS). For the opposite limit, that is a 1, w becomes flat and corresponds 
to the case of isotropic Large Angle Scattering (LAS). In §4.11 we will modify this functional 
form to introduce the possibility of anisotropic scattering. 



3. Deflection by a regular magnetic field in the upstream region 

It is well known that particle acceleration at a shock front with parallel magnetic field 
without scattering centers does not work. This magnetic scattering may be self-generated 
by the same particles, but the process of generation depends on the conditions in specific 
astrophysical environments. The case in which a regular magnetic field not parallel to the 
shock normal is present in the upstream fluid is quite interesting in that it allows for the 
return of the particles to the shock front even in the absence of scattering. In this section 
we investigate in detail the process of acceleration at shocks with arbitrary velocity when 
only a regular large scale magnetic field is present upstream (no scattering). We assume 
that enough turbulence is instead present in the downstream plasma to guarantee magnetic 
scattering of the particles. 

There are two main differences introduced by this situation when compared with the 
standard case considered in the previous section: 

• (a) Particle motion in the upstream region is deterministic: the stochasticity intro- 
duced by the interaction with scattering centers is assumed to be absent. This requires 
a new determination of the return probability P u introduced above. 

• (b) The presence of regular magnetic field with arbitrary orientation breaks the axial 
symmetry around the shock normal. This, in principle, would force us to treat the 
problem in the four angular variables jjL , <fio, /■* and (p. 



In the following we will show how addressing point (a) in fact solves point (b) as well. 
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3.1. Upstream return probability 

Let us Consider a particle entering upstream in the direction identified by the two angles 
fi and 0o, and returning to the shock along the direction identified by \i and 0. Since the 
motion of the particle is deterministic, the return direction is completely defined by the 
incoming coordinates, and we can write in full generality: 

Pu{fJ>o, 4>o] /i, 0) = {27fy 1 5 - /il(/i , <f>o)) 5 (0 - 0i (/i D , O )) , (11) 

where \i\ and 0i are obtained from the solution of the equation of motion, as discussed 
below. One can see that P u is effectively a function of only two variables. 

In order to apply the same mathematical procedure introduced in §2j we need to write 
P u as a function of azimuthal angles only. Therefore we use the properties of the delta 
function in S(/jl — /ii(/i , 0o)), to write: 



P u (/i o ,0 o ;/i,0) = — 
Ztt 



<90 o (/Z o ,/l) 



5 (0o - 0o) S (0 - 0i) 



dfx 

= P u (fi , n) 5 (0 O - O ) 5 (0 - 00 . (12) 

We now show that P tt (/i ,/i), as defined by Eq. ( [121) . is exactly the function to be used 
in eq. ([7]). This is easily shown by writing the fluxes of particles ingoing and outgoing the 
upstream plasma: 

J+{n,4>)= I dfi! d(f)' P u (/i',<f)',ii,<f))J-(n',(f)') , (13) 



J -i JO 
which, when integrated over 0, yields 



J+(/i) = / d0J+(/i,0)= / dfi' P u (//,//)J_(//) , (14) 



i 



where we as sumed that J_ is independent of 0. This is exactly the same relationship as 



was used in (jVietri 20031 ). and proves our point that the system may, in the average, still be 



treated as if it were symmetric about the shock normal. 

The key assumption here is that the flux crossing back into the upstream region from 
the downstream one, J_, be independent of the azimuthal angle 0. This is of course true in 
the Newtonian regime, because there the residence time for all particles diverges, and there 
is time for deflections to effectively erase anisotropies in the direction. But this must be 
true a fortiori in the relativistic regime, when one considers that the properties of scattering 
are of course still the same as in the Newtonian regime, while the surface to be recrossed, 
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i.e., the shock, is running away from the particles at a speed that becomes, asymptotically, 
a fair fraction of the particles' speed. So, while not exactly true, the independence of J_ 
from (f) is at least a good approximation. 

In order to write P u (/j, , fx) in a more explicit way, we need to solve the equation of 
motion of the particles, namely find the direction at which the particles re-cross the shock 
front as a function of the incoming direction. Particles move following a helicoidal trajectory 
around the magnetic field direction, indicated here as z. The problem is simplest if expressed 
in the frame O comoving with the upstream fluid but with the polar axis coincident with z. 
We mark with a tilde all quantities expressed in this frame. The equations of motion in the 
frame O are: 

jx{t) = jl , (15) 

4>(t) = 4> + ut, (16) 

where t is time and u is the Larmor frequency. The particles re-cross the shock when 
-ZparticicW = -ZshockW- This condition expressed in the frame O reads 

• / , , m .7 h cosa + (3 S 

sin(wi + (p ) — sin O = ~ — out, (17) 

sin a sin O 

where a is the angle between the shock normal z and the magnetic field direction z. The 
solution of Eq. ffTTj) gives the upstream residence time t* of the particles, to be evaluated 
numerically. 

The angles that identify the re-crossing direction, as functions of the residence time, are 

jxi = jl , and (18) 

01 = 4>o + Ut*(jjL ,(j) ) . (19) 

A rotation by the angle a provides us with the re-crossing coordinates fii and 0i in the fluid 
frame. At this point the Jacobian in Eq. (jT2l) can be calculated, although some care is needed 
because this Jacobian is not a single valued function: for each pair (/x ,aO the Jacobian has 
two values. This degeneracy arises because of the substitution of <p with /i, since each \x 
corresponds in general to two possible values of O - This is clear from fig. [Tj, where we show 
some examples of solutions: the directions of entrance and escape from the upstream fluid 
are plotted for different values of the shock speed and for different orientations of the large 
scale magnetic field. 

Eq. (fT7| admits a solution t* > only if the two following conditions are fulfilled: 



i) the initial velocity of a particle along the shock normal must be larger than the shock 
speed (otherwise the particle is prevented from crossing the shock to start with). This 
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Fig. 1. — Location of the particles entering the upstream region (dashed lines) and returning 
to the downstream region (solid lines) after being deflected by the magnetic field upstream. 
The directions are plotted in the plane T s f3 p ^ x - T s f3 P)y . (3 PtX and /3 Piy are the components 
of the particle velocity along the x and the y axis respectively. The origin corresponds to 
particles entering along the shock normal. Circles correspond to particles having constant 
fMo. The upper panels refer to T S /3 S = 0.04. The lower panel refer to T s j3 s = 1.0. In both cases 
we show the effects of three different orientations of the magnetic field (from left to right 
cos a = 0.0,0.4,0.8). The presence of an emp ty region when cos a > (3 S is due to particle 
leakage from upstream [Compare with fig. 1 in (lAchterberg et al. 20011 ) .] 



implies: 



jJ < -/?,, 



(20) 



ii) The particle velocity along the shock normal has to be less than the shock speed, 
namely 

jl cos a > —(3 S . (21) 

Particles not satisfying this last condition escape the shock region towards upstream 
infinity, a situation which is not realized in the case of scattering considered in £j2j This escape 
process occurs only for cos a > —/3 8 , and results in the loss of particles having the entrance 
pitch angles cosine exceeding /iminO^o, 4>o)- In fact for cos a < —(3 S , /z m i n = const = — 1 and 
all particles eventually re-cross the shock. 

When the particles are allowed to escape upstream, the acceleration is clearly expected 
to become less efficient and give rise to softer spectra of the accelerated particles (see §3.2p . 

Putting together all of the above, we can finally write the upstream conditional proba- 
bility as 

Ztc z — ' 

i=l,2 

where the sum is extended over the two branches of the Jacobian. 

For cos a < —f3 s the particles always return to the shock front and this forces the return 
probability to be unity when integrated over all outgoing directions: 

/ dnP v (n ,ii) = l. (23) 

J —u 

This integral condition is trivially satisfied by Eq. f)22p and is used as a check for P u after 
its numerical computation. 

Figs. El [3] and 0] show some examples of our calculations of P u (/x ,/i) as a function 
of \x for different values of /i D , for a Newtonian, a trans-relativistic and a relativistic shock 
respectively. For each case we show the results for different inclinations of the magnetic field 
with respect to the shock normal. It is worth noticing that P u does not change significantly 
when the inclination of the magnetic field varies in the range < cos a < /3 S , at a given 
shock speed. Therefore we do not expect a significant variation of the spectral slope in this 
range. In §3.21 we show that this is in fact the case. 



— - 8(-Ho -P s )Q{ll - fJtaOnijjLo, fi)) , (22) 
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Fig. 2. — Conditional probability P u (fi ,fi) as a function of the outgoing direction //, for a 
fixed value of the shock speed (T S (3 S = 0.04) with three different inclinations of the magnetic 
field (cos a = 0.0,0.4,0.9). For each plot the different lines correspond to different values of 
the ingoing direction fi . 




-1 -0.8 -0.6 -0.4 ^ -0.2 0.2 -1 -0.8 -0.6 -0.4 -0.2 0.2 -1 -0.8 -0.6 -0.4 ^ -0.2 0.2 




Fig. 4. — Like figj2]but for a relativistic shock (r s /3 s = 5.0; (3 S = 0.98). The three plots are 
very similar to each other because the condition cos a > f3 s is never reached. 
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3.2. Spectrum and Anisotropy of the accelerated particles for a large scale 

magnetic field upstream 

In this section we use Eq. and Eq. ([7]) to calculate the spectrum and angular 
distribution of the accelerated particles at the shock front. The return probabilities are 
calculated assuming that in the downstream fluid there is isotropic scattering, so that Pd 
can be calculated from Eq. (j3J) using Eq. (TTUT) as a scattering function. We assume a = 0.01 
for the SPAS regime and a = 10 for the LAS regime. In the upstream fluid we assume 
that particles can only be deflected by a large scale coherent magnetic field with arbitrary 
orientation with respect to the normal to the shock front. The return probability P u is 
therefore calculated as discussed in detail in §3.11 

The only information still lacking to proceed further is an equation of state for the 
medium, that would allow us to compute the vel ocity of the d ownstream fluid from the 



jump conditions at the shock front (see for instance (IGallant 20021 )). We assume that the gas 
upstream has zero pressure. Moreover, in the following we assume everywhere that the mag- 
netic field has no dynamical role, so that the standard jump conditions for an unmagnetized 
shock can be adopted (the role of the magnetic field becomes importa nt when the magnet ic 
energy density becomes comparable with the thermal energy density kirk fc Duffy 1999h ). 



Following much of the previ ous literature, we adopt the Synge equation of state for 



the downstream gas (ISynge 19571 ). assuming that only protons contribute. Although used 
widely, this assumption may not be well justified in a general case. We will illustrate our 
conclusions on the role of the equation of state for the spectrum and anisotropy of the 
accelerated particles in a separate paper. 



own 



Within this set of assumptions it is worth reminding that the compression ratio u up / u^, 
tends asymptotically to 4 for a non relativistic shock (even for shock speeds that are known 
to give lower compression factors) and to 3 for ultra-relativistic shocks. 

The simplest case to consider is that of a shock in which the large scale coherent mag- 
netic field in the upstream region is parallel to the shock front (cos a = 0). This is known as 
a perpendicular shock. The angular distribution and the slope of the spectrum of the accel- 
erated particles are plotted in Fig. [5] (the LAS (SPAS) case is shown in the left (right) panel) 
and Fig. O respectively, for various shock velocities ranging from newtonian to relativistic. 

The angular distribution of the particles in the downstream frame is seen to be rather 
anisotropic for the SPAS case, even in the newtonian regime. Large angle scattering (LAS) 
is evidently more efficient in isotropizing the accelerated particles. The anisotropies do not 
seem to affect the spectrum of the accelerated particles in the case of non relativistic shocks: 
the slope of the spectrum for both SPAS and LAS is 4.000 ± 0.001. The effect becomes 
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more prominent for faster shocks and in particular for relativistic shocks. In the SPAS case, 
for = 10, we found s = 4.272 ± 0.001, compatible with s = 4.28 ± 0.01, obtained by 



( lAchterberg et al. 200ll ) for T s = 10, with a Monte-Carlo simulation. 



In Fig. El the dotted and dashed lines refer to the SPAS and LAS cases respectively. At 
first sight it may appear rather surprising that in the limit of relativistic shocks the spectrum 
of accelerated particles is softer in the LAS regime than it is in the SPAS regime, since LAS 
is envisioned as more efficient in redirecting the particles to the shock front. This intuitive 
vision turns out to be incorrect, as also shown in Table 1, where we list the slope, the average 
energy gain and the return probability from downstream (as defined in Eqs. (1261) and fl29l) ) 
for a relativistic shock with T S /3 S = 5.0. 

One can see that while the average energy gain is similar in the two cases, the return 
probability in the case of LAS is 20% lower than for the SPAS case. Qualitatively this can 
be understood as follows: when the shock velocity increases, particles are caught up by the 
shock front when they have travelled only a small fraction (of order l/r s ) of their gyration. 
Once downstream, LAS is likely to swing them far from the shock front in a few interactions, 
while SPAS deflects their trajectories rather slowly yet remaining in the vicinity of the shock 
surface. This is responsible for the 20% difference in the average return probabilities in the 
two cases. This is also shown in Fig. H, where we plot the particles flux, J(/i) = \fi + Ud\g(fi) 
in terms of downstream coordinates: the total flux of particles entering the downstream 
section (— < \i < 1) is normalized to unity. It is clear from Fig. [7] that the flux of 
particles returning to the shock is slightly larger for the case of SPAS (dashed line in the 
range — 1 < ji < —Ud). 

A more interesting question concerns the effect of the orientation of the large scale 
magnetic field with respect to the normal to the shock. We have already emphasized that for 
any orientation different from that of a perpendicular shock, and in the absence of scattering 
processes upstream, particles are lost from the upstream region, because the shock cannot 
catch up with their motion. This happens when cos a > — /3 S , so that the phenomenon is 
increasingly more important for shocks approaching the parallel configuration. This reflects 



Table 1: Exact spectral slope, mean amplification and downstream return probability (as 
defined in eq. (12^|) and (1211 respectively) for T s j3 s = 5.0. 





slope 


(G) 


^p(downK 


SPAS 


4.218 ±0.001 


2.0387 


0.4165 


LAS 


4.445 ± 0.001 


2.0753 


0.3430 
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Fig. 5. — Particle distribution function at the shock front when a large scale coherent mag- 
netic field is present in the upstream region, with a direction parallel to the shock plane. In 
the downstream region particles are scattered in the LAS (left plot) or in the SPAS regime 
(right plot) (here the maximum of is arbitrarily set equal to 1). Several values of shock 
speeds are shown. The particle distribution functions always show a jump at /i = —(3 S . 
Large angle scattering makes distribution functions flatter compared with the small angle 
scattering case for — 1 < fi < —u d . 
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Fig. 6. — Spectral index vs. shock speed for the same configuration as in figJHJ 
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Fig. 7. — Particle flux J(/i) = \/i + -u rf |g(/z) across the shock front, as it appears in the 
downstream frame, when the upstream magnetic field is parallel to the shock (cos a = 0) 
and the downstream fluid is in the LAS regime (solid line) or SPAS regime (dashed line). 
The shock speed is T s j3 s = 5.0. The flux entering downstream (i.e. for — Ud < // < 1) is 
normalised to 1. In this way we see that downstream return probability, i.e. the integrated 
flux for — 1 < ji < —u d , is larger when the downstream region is in the SPAS regime. 
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in increasingly softer spectra. In the limit cos a — ► 1, all particles escape from the upstream 
region and no acceleration takes place. 

The slope of the spectrum as obtained from our calculations is plotted in Fig. [S] (solid 
lines and symbols) as a function of cos a for three different shock speeds (T S (3 S = 0.6, 1.0, 2.0): 
when there is no particle escape, the slope s is actually a constant, while it increases dramat- 
ically (and in fact diverges, showing the disappearance of the acceleration process) for values 
of cos a larger than —f3 s . In the small panel in Fig. [S]we also plot the return probability 
from upstream: for very inclined shocks the return probability is still very close to unity, as 
in the case of upstream scattering, but it drops rapidly for increasingly less inclined shocks. 

The steepening of the spectrum due to leak age of the particles towards upstream infinity 



can also be understood in terms of a Bell-like (IBell 19781 ) calculation, when carried out for 
the case of a large scale coherent magnetic field. The slope of the spectrum is related to the 
average return probability and to the average energy gain of the particles per cycle back and 
forth through the shock front through the expression: 

log(G) ' {M} 
where (G) is the mean amplification in a single cycle (downstream — > upstream — > down- 
stream), and (P re t) is the mean probability of returning to the shock. One should keep in 
mind that Bell's method, as expressed through the equation above is flawed in that it does 
not take into proper consideration the correlation between the amplification factor and the 
return probability. Moreover, Eq. (|24p hides the assumption of isotropy of the distribution 
function of the acce lerated particl es, since that formula was conceived in a discussion of non 



relativistic shocks ( (IPeacock 198ll ) introduced this formalism for particle acceleration at rel- 
ativistic shock fronts). All these limitations become of p articular imp ortance for relativistic 
shocks. A general expression for the slope was found in (jVietri 20031 ) . and reads: 

(P TCt )(G s - 3 ) = l. (25) 



In the following we use Eq. fl24j) . since we only want to provide the reader with an argument 
of plausibility for the steepening of the spectra in those cases in which particle leakage 
can take place in the upstream region. In order to account for this leakage, which cannot 
take place in the standard scenario of diffusive particle acceleration at a shock front, we 
generalize Eq. (12^1) in order to include the probability of escape from the acceleration box 
from upstream. This is easily achieved by replacing (P re t) with (P^) ■ (P^ ™ 11 ^). These 
mean values expressed in the downstream frame are: 

/p (down h _ /rr^° S- Ud d»9(n)(Ud + (J.)Pd(lJ.,lJ.o) 

Viet ) — r i , , 7z : 

j- Ud df*gw(u d + fi) 
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and 

/d(up)\ _ J- Ud d V III" dtiog(Vo)(u d + ti )P u (vo, fi) 

\ "ret / — r -u d , , w , s • \ Z ') 

J_! dn g(fi )(u d + Ho) 

In the last equation P u has also to be computed in terms of quantities evaluated in the 
downstream frame. Energy amplification for a particle entering the upstream region with 
direction /x (as measured downstream) and returning with direction fi, is obtained combining 
two Lorentz transformations: 

G(/i , /z) = 7 r 2 el (l - M re i/io) (1 + u ie ip) , (28) 

where /2 = (// + M rc i)/(1 + u re i/-t) is the returning direction as seen in the upstream frame. 
Averaging the amplification we have: 

_ /J" dfi g(fi )(u + /i ) a> G(no, /j)P u (/a , /x) 
/J" d/j, g(/j, )(u + (i a ) f^dfi P u (ji , fi) 

The spectral slope as computed through Eq. (12^1) is plotted in Fig. [8] (large box) with 
dashed lines; the corresponding upstream return probability (P^ ) is plotted in the small 
box (dashed lines). The agreement with our exact results is better than 1%, proving that the 
reason for the softening of the spectra of accelerated particles is in the increased probability 
that the particles leave the acceleration region when only a large scale coherent magnetic 
field is present upstream. 

The results discussed above apply to situations in which the magnetic field in the up- 
stream region can be considered as coherent on spatial scales exceeding the size of the accel- 
eration box. If the coherence scale of the field is smaller than the size of the accelerator, then 
the direction of the particles suffer a random wandering motion and one can think of this 
structured field as the source of diffusion and as a physical mechanism that imposes a max- 
imum energy to the accelerated particles (at least in the absence of radiative energy losses). 
Particles that escape from the shock region too fast (highest energy ones) have enough time 
to feel the effect of a coherent scale, while lower energy particles live in the accelerator 
for longer times and in principle may feel different orientations of the upstream magnetic 
field. This scenario is basically equivalent to having some degree of scattering upstream, an d 



should be treated with the formalism already discussed in (IVietri 20031 ; iBlasi fc Vietri 20051 ). 
As soon as a phenomenon equivalent to scattering is present, the probability of escape to 
upstream infinity vanishes, for all those particles that are confined in the accelerator for 
sufficiently long times. Moreover, one should keep in mind that even if a large scale coherent 
magnetic field is present to start with, the propagation of the accelerated (charged) particles 
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Fig. 8. — Large box: spectral slope as a function of cos a for three different values of the 
shock speed. The dashed lines show the spectral slope computed with Bell's method. Small 
box: the corresponding upstream return probabilities. 



in the upstream plasma is ver y likely to excite fluctuations in the magnetic field structure 
through streaming instability (jBell 19781 ). These fluctuations act as scattering centers and 
enhance the probability of returning to the shock front. 



4. Anisotropic scattering 



In this section we consider again the standard case in which particle motion in both 
the upstream and downstream fluids is diffusive, due to the presence of scattering agents. 
However, we include the possibility that the scattering, though spatially constant, may be 
anisotropic. The physical motivation for this generalization is the following: in a background 
of Alfven waves with a power spectrum Pw(k) (such that Pw{k)dk is the energy density in 
the form of waves with wavenumber in the range dk around k) the particles suffer angular 
diffusion with a diffusion coefficient 



where k r = fl/vfj, is the resonant wavenumber and Q is the gyration frequency of particles 
with momentum p in the background magnetic field Bq. One can clearly see from Eq. ( 1301) 
that the diffusion is anisotropic in general, unless the power spectrum has a specific ad hoc 
form. One should keep in mind that Eq. (!30l) is obtained in the context of quasi-linear 
theory. A full non-linear treatment might show how the turbulence is distributed and which 
is the resulting particle angular distribution. 

In the calculations that follow, we quantify the effects of anisotropic scattering on the 
spectrum and angular distribution of the accelerated particles. The calculation of specific 
patterns of anisotropy in the scattering agents is beyond the scopes of this paper, therefore 
we adopt a few simple but physically meaningful toy models of anisotropic scattering and 
we carry out the calculations within those models. 



We parametrize the anisotropy in such a way to reproduce the following four patterns: 




(30) 



4.1. Modelling anisotropy 



case A: Particles are scattered per unit length more efficiently while they move away 
from the shock front than they are on their way to the shock front, both upstream and 
downstream. 
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• case B (opposite of case A): Particles are scattered per unit length more efficiently on 
their way to the shock front than they are while they move away from the shock front, 
both upstream and downstream. 

• case C: In the downstream fluid, particles are scattered per unit length more efficiently 
while they move away from the shock front (p —*■ 1) than they are on their way to the 
shock front (//—>• —1). In the upstream fluid the situation is reversed, and scattering 
is more efficient for the particles that are moving toward of the shock [jj, — > 1). 

• case D (opposite of C): Scattering is more effective around \i ~ — 1 both upstream 
and downstream. 

A pictorial representation of cases A-D is shown in Fig. [91 

In order to simulate the cases A-D above, we adopt a scattering function similar to Eq. 
(ITO]) . but modified to introduce anisotropic scattering. In particular, to achieve this goal we 
allow the width a of the scattering function to depend on both the initial and final directions 
/jf and [A, so that: 



It is worth stressing that the scattering function has to be symmetric if we exchange \i 
with // as a consequence of Liouville's theorem, so we are forced to look for a symmetric 
function cr(yU, /i'). 

In order to apply the functional form Eq. (131]) to the cases A-D, it is sufficient to adopt 
the following expression for the mean scattering angle a(/J,,f/): 



Both cr_|_ and cr_ have o"o as the maximum and cr /a as the minimum value. For this reason 
we will refer to a as the Anisotropy Factor. For a = 1, isotropic scattering is recovered. 

The resulting scattering function if T (/i, //), obtained substituting Eq. (|32|) into (I3TI) . is 
plotted in Fig. [TU1 together with the isotropic scattering function (Eq. ({TO]) ), for a = 0.05 
and a = 10. These plots clarify how w + and u>_ can simulate a scattering more efficient in 
the fi = +1 and /i = — 1 directions respectively. 

The condition f_ t dfiw(fi, //) = 1 that states the probability conservation, is fulfilled by 
Eq. fl3TT) provided o"o 1. In the numerical calculations that follow we assume o"o = 0.05. 




(31) 




(32) 
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Fig. 9. — Pictorial representation of the four patterns of anisotropic scattering considered in 
our calculations. 
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Fig. 10. — The thick lines show the anisotropic scattering functions /i ) (left box) and 

W-(fi,fi ) (right box), as functions of /i and for different values of the incoming direction 
/i . The anisotropy factor is a = 10 and a = 0.05. For comparison the isotropic scattering 
function (Eq. ( flOi) ) is shown with thin lines and for the same values of /i . 
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Using a- and a + in different combinations for the upstream and the downstream fluids, 
we can reproduce scenarios A, B, C, and D, as summarized in Table [2j 



4.2. Results: anisotropic scattering for shocks of arbitrary speed and fixed 

anisotropy factor 

Following the procedure outlined in §|2]and making use of Eqs. (|3"Ij) and (132]) . we compute 
the spectral index and the angular distribution for the scenarios A, B, C and D, described 
above. In each case both the parameter <j and the anisotropy factor a are fixed (a = 0.05 
and a = 10), while the shock velocity is allowed to vary within the range 0.04 $C T S /3 S ^ 5. 

The angular part of the distribution function is shown in Fig. [11] for the scenarios A 
(left panel) and B (right panel) and in Fig. [12] for the scenarios C (left panel) and D (right 
panel). The slope of the spectrum of accelerated particles is plotted in Fig. [T31 

For relativistic shocks, the spread in the slope of the spectrum of accelerated particles 
has less spread around ~ 4, although in general it remains true that harder spectra are 
obtained in the scenarios B and D. 

A note of caution is necessary to interpret the apparent peak in the slopes at T S /3 S ~ 1 
for cases A, and at T S (3 S ~ 3 for cases D. These peaks are completely unrelated to anisotropic 
scattering and is instead the result of the bre aking of the regime of small pitch angle scattering 



or SPAS), as was already pointed out in (IBlasi fc Vietri 20051 ). The acceleration process 



does no longer take place in the SPAS regime when > l/4er, which happens at higher 
Lorentz factor when a is smaller. This is shown in Fig. [T4J where we plot the slope for 
the case of isotropic SPAS for a = 0.1 (dashed line) and o = 0.01 (solid line), and th e 



corresponding angular distribution for T S (3 S = 5. As already found in (IBlasi &: Vietri 20051 ). 
the transition from SPAS to LAS is generally accompanied by a hardening of the spectra 
of accelerated particles. The peak seen in Fig. [T3] is simply the consequence of an effective 
value of a in the anisotropic scattering cases A and D. This is also clear comparing angular 
distributions of Fig. [TH with the angular distribution of cases A and D for T S (3 S = 4 and 
5: the curves show the same behaviour with a jump at \i — — /3d OW n and a peak that moves 

Table 2: Summary of mean scattering angle used in the different scenario of figM 
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Fig. 11. — Particle distribution function at the shock front for anisotropic scattering of type 
A (on the left) and B (on the right) both with a = 10. Each line represents a different shock 
speed as the labels show. 



Fig. 12. — The same as figOJ] but for anisotropic scattering of type C (on the left) and D 
(on the right). 
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towards jj, — 1 as the shock speed increase. 



5. Conclusions and Discussion 

In this paper we carried out exact calculations of the angular distribution function and 
spectral slope of the particles accelerated at plane shock fr onts moving with arbitrary velocity , 



generalizing a method previously described in detail in (jVietri 20031 ; iBlasi &: Vietri 20051 ). 
In particular, we specialized our calculations to two situations: 1) presence of a large scale 
coherent magnetic field of arbitrary orientation with respect to the shock normal, in the 
upstream fluid; 2) possibility of anisotropic scattering in the upstream and downstream 
plasmas. 

Our calculations allowed us to describe the importance of the inclination of the magnetic 
field when this has a large coherence length and there are no scattering agents upstream. 
For newtonian shocks, only quasi-perpendicular fields (namely perpendicular to the shock 
normal) are of practical importance, in that the return of particles to the shock from the 
upstream section is warranted. Quasi-parallel shocks imply a very low probability of return, 
so that the spectrum of accelerated particles is extremely soft. The process of acceleration 
eventually shuts off for parallel shocks. For relativistic shocks, the situation is less pessimistic 
because the accelerated particles and the shock front move with comparable velocities in the 
upstream frame. In general, the acceleration stops being efficient when the cosine of the 
inclination angle a of the magnetic field with respect to the shock normal is comparable 
with the shock speed in units of the speed of light. The slope of the spectrum of accelerated 
particles for cos a = as a function of the shock velocity is plotted in Fig. for the two 
cases in which SPAS or LAS is operating in the downstream plasma. The slope as a function 
of cos a = for shocks moving at different speeds is shown in Fig. [HJ In the same figure 
we also show the return probability from the upstream section, in order to emphasize that 
the presence of a large scale magnetic field upstream leads to particle leakage to upstream 
infinity. This latter phenomenon disappears when scattering is present, in that scattering 
always allows for the shock to reach the accelerated particles. In this case the probability 
of returning to the shock at an arbitrary direction is unity. One can ask when and how the 
transition from a situation in which there is no scattering to one in which scattering is at work 
takes place. When some scattering is present but the energy density in the scattering agents 
(e.g. Alfven waves) is very low compared with the energy density in the background magnetic 
field, only very low energy particles are effectively scattered. When their energy becomes 
large enough, they only feel the presence of the coherent field. Increasing the amount of 
scattering, this transition energy becomes gradually higher. Particles whose Larmor radius 
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Fig. 13. — Slope vs. shock speed in the four different anisotropic scattering scenarios: A 
and B in the left box, C and D in the right one. All plots are obtained with a = 10 and 
(j = 0.05. Slope resulting from isotropic scattering (computed with a = 0.01) is also shown 
for comparison (solid line). 



4.2 




Fig. 14. — These plots show the breaking of the SPAS approximation. Left box: slope of 
the spectrum of accelerated particles in the case of isotropic scattering with o = 0.01 (solid 
line) and a = 0.1 (dashed line). Right box: angular distribution for T S (3 S = 5 with a = 0.01 
(solid line) and a = 0.1 (dashed line). 
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is larger that the coherence scale of the magnetic field can eventually escape the accelerator. 
In general the level of turbulence (and therefore of scattering) and the number of accelerated 
particles are not independe nt since the turbulence may be self-generated through streaming- 



like instabilities flBell 19781 ) 



In §4] we extended our analysis to the very interesting case of anisotropic scattering 
in both the upstream (unshocked) and downstream (shocked) medium. The pattern of 
anisotropy, which clearly depends on the details of the formation and development of the 
scattering centers, has been parametrized in four different scenarios, and for each one we 
calculated the angular part of the distribution function and the spectrum of the accelerated 
particles. Deviations from the predictions obtained in the context of isotropic SPAS and 
LAS have been quantified: the typical magnitude of these deflections is a few percent, but 
there are situations in which the deviation is more interesting, in particular because it goes 
in the direction of making spectra harder. 

This research was partially funded through grant COFIN 2004-2005. PB is grateful to 
the KIPAC at SLAC/Stanford for hospitality during the final stages of preparation of the 
manuscript. 
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